Annealing mechanisms of intrinsic defects in 3C-SiC: a theoretical study 
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The annealing kinetics of mobile intrinsic defects is investigated by an ab initio method based 
on density functional theory. The interstitial-vacancy recombination, the diffusion of vacancies and 
interstitials to defect sinks (e.g. surfaces or dislocations) as well as the formation of interstitial- 
clusters are considered. The calculated migration and reaction barriers suggest a hierarchical or- 
dering of competing annealing mechanisms. The higher mobility of carbon and silicon interstitials 
as compared to the vacancies drives the annealing mechanims at lower temperatures including the 
vacancy-interstitial recombination and the formation of interstitial carbon clusters. These clusters 
act as a source for carbon interstials at elevated temperatures. In p-type material we discuss the 
transformation of the silicon vacancy into the more stable vacancy-antisite complex as an annealing 
mechanism, which is activated before the vacancy migration. Recent annealing studies of vacancy- 
related centers in irradiated 3C- and 4H-SiC and semi-insulating 4H-SiC are interpreted in terms of 
the proposed hierarchy of annealing mechanisms. 

PACS numbers: 61.72. Ji,61.72.Cc,66.30.-h 
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I. INTRODUCTION 

The unique features of silicon carbide, such as the wide 
band gap and the electrical and thermal stability, recom- 
mend this semiconductor for high power, high frequency 
and high temperature applications. As for any semicon- 
ductor, unwanted defects are introduced into this mate- 
rial by growth processes and ion- implantation of dopants. 
Besides the primary defects - interstitials and vacancies 
- secondary defects such as antisites and defect clusters 
are created. Thermal annealing is applied to reduce these 
defects to an inevitable abundance. Thereby the high 
electron mobility is restored and the electrical activation 
of dopants is achieved. 

The annealing properties of intrinsic and impurity- 
related defect centers were studied in irradiated SiC by 
various experimental techniques, for instance electron 
spin resonance techniques (EPR), deep level transient 
spectroscopy (DLTS), photoluminescence spectroscopy 
(PL) and positron annihilation spectroscopy (PAS). 
Some of the reported centers persist up to a tempera- 
ture range of 1300°C-1700°C. Examples for such defects 
arc the PL-centers Di (Refs. IHIH and D„ (Refs. |H 
as well as the DLTS-center Z 1 /Z 2 and Ei/E 2 in 4H- 
and 6H-SiC (Ref. 0) respectively. The microscopic ori- 
gin of most of the centers is still investigated. The PL- 
centers P— U (Ref. 0) and Dn have been interpreted as 
carbon interstitial^ or interstitial clusters^ due to their 
carbon-related localized vibrational modes (LVMs) with 
frequencies above the SiC bulk-spectrum. Depending on 
the experimental conditions some of these centers even 
grow in intensity during the heat treatment in a tem- 
perature range between 500°C and 1200°C, when others 
(e.g. vacancy-related centers) gradually vanishA^ 

Vacancies and interstitials act as vehicles for other- 
wise immobile point defects. Defect clusters may only 
be diminished by reacting with or separating into these 
mobile entities. Vacancy related defects have been stud- 



ied by PAS 6 ' 10 ! 11 ! 12 ! 13 and EPR 14 i 15 i 16 i 17 in irradiated 
material. In PASfi - defects related to silicon and car- 
bon vacancies can be distinguished thanks to a consid- 
erable difference of the positron lifetimes, as predicted 
by theoryii2*iS EPR-centers have also been identified as 
isolated siliconi2*i£*i£ and carboniSiSSiSi vacancies and 
the assignment has been verified theoreticallyii^i 2 ^* 2 ^ 2 ^ 
Characteristic differences were observed for the anneal- 
ing behavior of carbon and silicon vacancies in irradi- 
ated SiC. A lower annealing temperature (^500°C) was 
deduced for the carbon vacancy by PAS 10 and EPR. 16 
The annealing of the silicon vacancy occurrs in several 
stages, at 150°C, 350°C and 750°C (the final anneal- 
ing stage) as observed by EPRii in 3C-SiC in agree- 
ment with PASi^ data. In irradiated n-type 4H- and 
6H-SiC, however, additional annealing stages of centers 
related to a silicon vacancy were observed by PAS&2^ at 
temperatures up to 1450°C. This finding was explained 
by the formation of stable nitrogen-vacancy complexes. 
On the other hand, a vacancy-antisite complex was iden- 
tified by EPR experiment* 26 in irradiated n-type 6H- 
SiC annealed above 750° C. An EPR-center with simi- 
lar properties was also identified in 3C-SiCi^£ The cen- 
ter evolves from the silicon vacancy as a consequence of 
the vacancy's metastability in p-type and intrinsic ma- 
terial predicted previously by theory^Si In both exper- 
iments the signature of the silicon vacancy was absent. 
Other experiments 20 ^ indicate that the carbon and sil- 
icon vacancies posses a higher thermal stability in semi- 
insulating than in irradiated SiC. A comprehensive inter- 
pretation of the annealing experiments in terms of simple 
mechanisms is lacking sofar. In particular, the role that 
interstitials and vacancies play in the reported annealing 
stages is not known. Furthermore, the Fermi-level effect 
on the annealing behavior of the charged defects has not 
been clarified. 

In the present paper we analyze the annealing of va- 
cancies and interstitials in 3C-SiC at a microscopic level. 
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We study the effect of clustering on the defect kinetics in 
case of carbon interstitials. An ab initio method within 
the framework of density functional theory (DFT) is em- 
ployed to calculate the migration barriers and recombi- 
nation paths. The aggregation of carbon intcrstititials is 
analysed and the dissociation energy needed to re-emit a 
single carbon interstitial is calculated. A hierarchy of an- 
nealing mechanisms (cf. Fig. ^| is proposed on the basis 
of the calculated barriers and dissociation energies at dif- 
ferent doping conditions. For the carbon-related defects 
we show that the vacancy-interstitial recombination and 
the diffusion-limited annealing of carbon interstitials are 
the first two annealing stages in this hierarchy followed 
by the diffusion-based annealing of carbon vacancies. Be- 
sides the diffusion-limited vacancy-interstitial recombi- 
nation the clustering of carbon interstitials may facili- 
tate their annealing. Thermally stable carbon clusters 
re-emit carbon interstitials at high temperatures thereby 
contributing to the kinetics of other thermally stable de- 
fects. 

While the Fermi-level effect is not pronounced in the 
case of carbon vacancies and interstitials, it substantially 
affects the annealing of silicon vacancies and intersti- 
tials. This is related to the metastability of the silicon 
vacancy in p-type and intrinsic material and the exis- 
tence of two different interstitial configurations in p-type 
and intrinsic or n-type SiC. For intrinsic (compensated) 
SiC, the vacancy-interstitial recombination preceeds the 
metastability-related transformation of the vacancy into 
the carbon vacancy-antisite complex, which in a next 
stage either dissociates or anneals via a diffusion-limited 
mechanism. For p-type conditions the transformation- 
induced annealing dominates over the diffusion-limited 
vacancy-interstitial recombination. In n-type material a 
diffusion-limited annealing of the silicon vacancy should 
follow the vacancy-interstitial recombination in the hier- 
archy of annealing mechanisms. 

The outline of the paper is as follows. Preceeded by 
a description of the method in Section [H] we summarize 
the migration mechanisms of interstitials and vacancies 
in Section IIIII We also review the metastability of the 
silicon vacancy there and analyze the kinetic aspects of 
this metastability with emphasis on Fermi-level effects. 
Section I1VI treats the stable Frenkel pairs and the re- 
combination paths of vacancies and interstitials. In Sec- 
tion the properties of carbon interstitial clusters and 
their dissoziation energies are described. The hierarchy 
of annealing mechanisms is outlined and discussed in the 
light of recent experiments in Section IVII 
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FIG. 1: Possible annealing paths of vacancies: (a) carbon va- 
cancy: recombination of a carbon split-interstitial (C ap / 10 o)) 
with a carbon vacancy (Vc) and the migration of carbon va- 
cancies to sinks, (b) silicon vacancy: Frenkel pair recombi- 
nation (compensated and n-type material) and diffusion lim- 
ited recombination with split interstitial (Si sp (no)), and (c) 
annealing by the metastability induced transformation into 
a carbon vacancy-antisite complex (p-type and compensated 
material) and the migration of silicon vacancies to sinks (n- 
type). 



their environment are described using large supercells, 
equivalent to 64 and 216 crystal lattice sites. A spe- 
cial k-point mesh*» with 8 k-points in the Brillouin zone 
(2x2x2-mesh) is used. For the description of Jahn- 
Teller distortions we performed the calculations in 216 
atoms cell at the T-point in order to maintain the de- 
generacy of the defect levels. The ionization levels of 
charged defects were calculated following the approxi- 
mate procedure proposed by Makov and Payne*2 to ac- 
count for the error introduced by the compensating back- 
ground. Especially for the highly charged interstitial de- 
fects this procedure leads to more consistent results iKBiiSi 
Norm-conserving carbon and silicon pseudopotentials of 
the Troullier-Martins typ o 40 ' 41 are employed. The car- 
bon pseudopotential has been optimized 42 for calcula- 
tions with a small basis set. Extensive tests showed that 
a basis set of plane-waves with a kinetic energy up to 
30 Ry yielded practically converged total energy differ- 
ences. The coordinates of all atoms in the simulation 
cell have been relaxed. Test have shown that all relevant 
relaxations are well contained within the simulation cells. 



II. METHOD 



A. Formation energies 



We employ the plane wave pseudopotential program 
package FHI96MD 31 > based on density functional the- 



,32.33 



and the local density approximations 4 ** 5 . (LDA) 
for the exchange-correlation functional. Spin effects are 
included within the LSDA where noted. Defects and 



The abundance of defects is determined by their ef- 
fective formation energy Ef. For charged defects, the 
latter depends on the doping conditions via the Fermi 
level fjL-p and on the stoichiometry of the material via 
the carbon and silicon chemical potentials /ic and /xsi- 
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In equilibrium nc and /xsi are not independent but re- 
lated to the chemical potential of SiC /Ltgic through 
MSiC = + A*Si- Choosing /isi as a free parameter, we 
express it as usi = Msi ^A* m terms of the chemical 
potential of crystalline Si and the difference A/i to this 
value. The foramtion energy is given by 

E{ = Er> t ce n- nsi^si—nc fJ-c — n e (Ev + fi F ) , (1) 

where -Ed, ceil is the total energy of a super cell containing 
the defect and up is defined relative to the valence band 
edge Ey< As outlined in Ref. |3^ Ef may be conviniently 
written as 

E { = Eg c - (n S i - n c ) - n c /i F . (2) 

In the above equations n c is the number of excess elec- 
trons at the defect, nsi and nc are the numbers of silicon 
and carbon atoms in the supercell. The quantity E^J is 
determined by E^ c = E^ ccll - n c MSiC - ("Si - nc) A*si _ 
n c Ey where the chemical potentials /isic and /Xgj are to 
a good approximation given by total energies of the per- 
fect crystals. For Ey we take the value from a defect- 
free cell including a correction for the difference of the 
average potential in the defect and defect-free cell (c.f 
Ref. 1391 . To ensure the stability of SiC, A^ may only 
vary between the heat of formation of the correspond- 
ing polytype Hf, gic and 0, reflecting C-rich and Si-rich 
conditions respectively. -fff,siC only slightly depends on 
the polytype and amounts to 0.61 eV according to our 
findings choosing diamond as a reference for uc . This is 
slightly lower than the experimental value of 0.72 eV. 

The ionization level of charged defects corresponds to 
the value of n-p at which the defect alters its charge state. 
Thus it determines the defect charge state for given dop- 
ing conditions. The ionization level is obtained from the 
formation energy by 

e(q 2 \ qi ) = Eg - Eg . (3) 

Here qi and q^ indicate the different charge states of the 
defect. Usually the charge state changes by one electron. 
However, in some cases^ it may change by two when 
the additional electron-electron repulsion can be compen- 
sated by a lattice relaxation. This effect, known as the 
negative-U effect, leads to an attractive effective electron- 
electron interaction so that the ionization level e(q+ l\q) 
appears below e(q\q— 1). The realization of a charge state 
as a ground state for implies that the corresponding ion- 
ization levels and the occupied single particle levels are 
located below the conduction band edge. Even though a 
physical relevance of the Kohn-Sham levels is not rigor- 
ously founded by DFT, they reproduce the experimental 
quasi-particle band structure quite well. Yet, the calcu- 
lated Kohn-Sham band gaps (1.2 eV and 2.2 eV in 3C- 
and 4H-SiC) are typically smaller than the experimental 
values (2.39 eV and 3.27 eV). In order to assess the defect 
charge states we follow a common practice and use the 
experimental value for the conduction band edge. How- 
ever this procedure cannot unambigiously describe the 



charge states that are stable only in the vicinity of the 
conduction band edge. 

B. Defect migration 

The analysis of migration and recombination paths 
starts with establishing the stable initial and final con- 
figuration. Then standard search methods are used to 
find the path between these configurations and to deter- 
mine the saddle point. For the interstitial migration and 
the Frenkel-pair recombination we employ an implemen- 
tation 44 of the ridge method of Ionova and Carter. 4 '"' The 
ridge method performs an automatic search for the sad- 
dle point starting from the initial and final configuration. 
We repeat the search for all relevant charge states and 
analyze the obtained saddle points with respect to the ge- 
ometry and the electronic structure. When the migration 
path depends on the occupation of defect states (near 
degeneracies, Jahn- Teller distorted ground states) this 
automatic search fails. This is relevant for the vacancy 
migration and the metastability-induced transformation 
of the silicon vacancy into a vacancy-antisite complex 
as discussed in Ref. yij In these cases, we analyze the 
potential energy surfaces using the drag method. The 
procedure is outlined in Ref. OS where we discuss the 
migration of the vacancies and their transition states in 
detail. 



III. MIGRATION OF INTERSTITIALS AND 
VACANCIES 

A. interstitials 

In 3C-SiC a large number of interstitial configurations 
exists. Besides the tetrahedral and hexagonal intersti- 
tials also split-interstitials are relevant. In the latter the 
interstitial atom and the lattice atom share a lattice site 
in a dumbbell-like configuration, which is preferentially 
oriented along the (100)- or (HO)-direction. Recently, we 
have analyzed these configurations for silicon and carbon 
interstitials^24& In particular, we have discerned the en- 
ergetically most favorable sites. At these sites the in- 
terstitial migration starts and ends. Other interstitial 
configurations act as intermediate states along the mi- 
gration path. The migration path and the energy bar- 
riers were investigated using a 64-atom cell and special 
k-points (2x2x2 Monkhorst-Pack mesh). 

Among the carbon interstitials the two split- 
interstitials C sp (ioo) and C sp si(ioo> are the most favorable 
configurations. The split-interstitial C sp ^ 10 o) is centered 
at a carbon site and oriented in the (lOO)-direction as dis- 
played in Fig.|2k. It is more favorable than the C sp si(ioo>- 
interstitial which is a (lOO)-oriented silicon-carbon dumb- 
bell at a silicon site. The electronic structure of the in- 
terstitials is discussed in detail in Ref. y!J The two inter- 
stitials are stable in positive and negative charge states. 
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TABLE I: Ionization levels of the mobile defects in 3C-SiC TABLE III: Migration of silicon interstitials: energy barriers 
given relative to the valence band edge . Em in eV for different migration paths. The dominant migra- 

tion paths are shown in Fig. [5] The charge state is indicated. 
See text for details. 
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e(l-|0) 


e(2-|l-) 
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1.29 


1.14 


2.69 


2.04 


V Si 




0.18 


0.61 


1.76 


Vc-Csi 


1.24 


1.79 


2.19 




Si sp {iio> 


0.4 


1.1 






C S p(100) 


0.6 


0.8 


1.8 




C sp si(ioo> 


0.4 


0.7 


1.9 


2.3 



51 TC^ O1 sp(100) O1 TC^ O1 sp(110) S1 TC ^ S1 TSi 

1- »2 3.5 3.4 3.5 

2- sd 0.03 0.05 0.3 

q-+ n-o a-o 

B1 sp{H0) O1 sp(110) al sp{110) 

Sl Jp(110) ^ Sl sp(110) ^ Sl sp{110) 



l<->2 1.0 1.4 1.3 



TABLE II: Migration of carbon interstitials: energy barriers 

of nearest neighbor hops and second neighbor hops for the 

~<o 

^sp{100) 



relevant charge states. The migration of C° /ln n\ starts from 



the tilted configuration. 



path 




energy barrier 


[eV] 






2+ 


1+ 





1" 


Csp(ioo) - * C sp si(ioo> 


1.7 


0.9 


0.5 


0.7 


C S pSi(100) — * C sp (ioo) 


0.9 


0.2 


0.2 


0.1 


C sp (ioo> C sp (ioo> 


1.4 


0.9 


0.5 


0.6 



The ionization levels of C sp ( 10 o) and C sp si(ioo) are given 
in Tab. [I] Other interstitial sites, such as the hexagonal 
site and the tetrahedral sites, are higher in energy. 

Two alternative migration paths have been found for 
the migration of carbon interstitials: a migration by sec- 
ond neighbor hops of C sp (ioo) between adjacent carbon- 
sites and a migration by nearest neighbor hops between 
adjacent carbon and silicon lattice sites, going through 



the sequence C sp ( 



100) 



pSi(100) 



p(ioo) 



The migra- 



tion by second neighbor hops is indicated in Fig. [5] 

The migration barriers for the two mechanisms 
are given in Tab. [n] The nearest neighbor hop 



C. 



p(ioo) 



C, 



pSi(100) 



and the second neighbor hop 



Csp(ioo) ^C sp (ioo) have similar migration barriers, ex- 
cept for the C sr J", 100 > . Small barriers are found for the mi- 
gration of the positive, neutral and negative interstitials. 
The analysis of migration paths via other intermediate 
interstitial sites indicates that such paths are unlikely: 
the energy differences between C sp (i o) and these sites 
typically exceed the small migration barriers found for 
the above mechanisms. 

Among the silicon interstitials the important sites are 
the tetrahedrally carbon coordinated site (Sixc) an d the 
(llO)-oriented split-interstitial (Si sp mo))- The tetrahe- 
dral interstitial dominates in p-type material. Our anal- 
ysis of the electronic structure shows 59 that it has no 
deep defect states within the band gap and therefore 
has a charge state of 4 + . The bare silicon ion is effi- 
ciently screened by a crystal polarization. Considering 
the possible occupation of a defect state located slightly 
above the experimental conduction band edge we con- 
cluded in Ref |3^ that the charge state 4 + prevails for 



Fermi-levels below the mid-gap (by the large Madelung 
corrections e(3 + |4 + ) is placed above the mid-gap). This 
is also consistent with our recent findings for Sixc in 4H- 



SiC. The split-interstitial Si 



sp(no) 



dominates in intrinsic 



and n-type material. It has several deep states within 
the band gap and is stable in the charge states 2 + , 1 + , 
and the neutral state. The ionization levels are given in 
Tab. [3 For the diffusion in intrinsic and n-type material 
the neutral configuration is most important. 

In p-type material all migration paths begin and end at 
the Sixc site. The migration of Sixc proceeds basically 
by two different mechanisms: (i) a kick-out mechanism 



via via the split-interstitials Si. 



sp (100) 



and Si 



.p<110> 



as in- 
termediate sites and (ii) a direct migration via a silicon 
coordinated interstitial site Sixsi- I n the kick-out mech- 
anism the tetrahedral interstitial moves towards one of 
its silicon neighbors and forms a split-interstitial as an 
intermediate site as indicated in Fig. [21 To complete the 
migration, one of the silicon atoms of the split-interstitial 
moves to a neighboring Sixc _ site. The energy barriers 
for the migration via the different intermediate sites are 
listed in Tab. IIIII The lowest barrier is obtained for the 
migration via Si sp /no\- The barriers for the other paths 
are only slightly higher 0.1 eV. In case of the kick-out 
mechanism the transition state lies in the vicinity of the 
intermediate site, which is found to be hardly stable. For 
the direct path the intermediate Sixsi is more stable with 
a barrier of 0.3 eV. 

In intrinsic (compensated) and n-type material the in- 
terstitial migration of Si° p ^ 110 ^ dominates. As indicated 
in Fig. |2t the migration proceeds by hops of one of the 
silicon atoms to the next silicon lattice site, where a new 
Si sp (no)-interstitial is formed. A barrier of 1.4 eV is found 
for this hop. The interstitial may also change its orienta- 
tion. For example, by the rotation of the pair around the 
(OOl)-axis the orientation changes from (110) to (110). 
This process is associated with a barrier of 1.3 eV. We 
have also considered the positive interstitial, which mi- 
grates with a barrier of 1.0 eV. 

So far the migration of interstitials has not been inves- 
tigated on an ab initio level in 4H-SiC. However, our anal- 
ysis of the ground state configurations shows^i that the 
carbon split-interstitials at the carbon and silicon lattice 



C sp <100> 



Si sp <no> 
Si TC Si TC 



Si S p<110> 



FIG. 2: Migration of interstitial: (a) migration of the carbon split-interstitial by second neighbor hops, (b) kick-out migration 
of the silicon interstitial between adjacent Si-rc-sites via the split-interstitial Si sp (no) , and (c) migration of the silicon interstitial 
between adjacent Si sp (no configurations. 



sites have identical properties as in 3C-SiC, except that 
^Wioo) becomes stable in n-type 4H-SiC. For the sili- 
con interstitials in 4H-SiC, besides the sites discussed for 
3C-SiC, additional sites appear: a site where the inter- 
stitial has simultaneously four carbon and silicon neigh- 
bors, and an open cage formed by the hexagonal rings 
of two adjacent lattice planes. For the interstitial sites 
with cubic environment we have found similar properties 
to 3C-SiC, also with respect to the (llO)-oriented split- 
interstitial. In the open cage a position close to three 
carbon neighbors is preferred. The interstitial with four 
silicon and carbon neighbors is energetically unfavorable 
and may represent an intermediate site of a migration 
path. From the bonding properties of the interstitials we 
expect similar migration barriers as in 3C-SiC. 



B. Vacancy migration 

The ground state properties of the carbon and silicon 
vacancy have been discussed in detail by Zywietz et til J* 
and Torpo et al. 49 . Here we briefly summarize a few 
facts for the sake of clarity. Both vacancies have a three- 
fold degenerate level within the band gap. Its occupation 
leads in case of the carbon vacancy (Vc) to a consider- 
able Jahn- Teller-relaxation. The vacancy can exist in the 
charge states 2 + , 1 + , and neutral. Our calculation indi- 
cate that the negative charge states should not be real- 
ized. TheionizationlevelsarelistedinTab.nl According 
to our calculations the positive charge state is unstable 
due to a negative-U behavior. However, for U , which is 
given by U — e(+\ + +) — e(0|+), we obtain a rather 
small value of — 0.15eV. This value should be considered 
as an estimate that depends on the size of the supercell 
and on the correction for the total energy of the charge 
defect^ Neglecting this correction, Zywietz et alM- ar- 
rived at U = — 0.48eV. In contrast to our findings, the 
calculations of Torpo et al~ suggest the absence of a 
negative-U behavior. However, the deviations in the re- 
sults for U are within the uncertainty in the evaluation 



of total energies for charged defects of about 0.2 eV. Irre- 
spective whether a small negative-U is present or not, the 
positive charge state becomes relevant for /ip around the 
ionization level e(0| + +) = 1.2 eV, in particular, when 
the Fermi level is pinned by an excess concentration of 
vacancies. This is also suggested by the recent identifi- 
cation of the EI5-center in 4H-SiC as a positive carbon 
vacancy^&2£ which shows a similar behavior as in 
3C-S1C. 

The DFT-LDA calculations predict that the silicon 
vacancy prefers a high-spin state over a Jahn- Teller- 
distorted ground state. This is consistent with the find- 
ings of spin-resonance experiments 1 ^*^ for Vg ; . Accord- 
ing to our analysis the silicon vacancy is realized in the 
charge states 1 + , neutral, 1~, and 2~ in 3C-SiC. Al- 
ready in p-type material the vacancy becomes neutral 
and changes to a negative charge state for a Fermi-level 
at 0.6 eV. The charge state 2~ is stabilized in n-type ma- 
terial. 

In 4H-SiC the carbon and silicon vacancy possess sim- 
ilar properties. The occupation of the lattice sites with 
local cubic and hexagonal coordination introduces neg- 
ligible changes for the silicon vacancy. In case of Vgj 
the calculation of HF-tensors 24 and EPR-experimentsi» 
show that the two sites cannot be distinguished. Also for 
the carbon vacancy similar formation energies and ion- 
ization levels are obtained at the cubic and hexagonal 
site. However, the different rehybridisation of the silicon 
dangling bonds 24 affects the local symmetry and for Vj 
leads to different HF-signatures at the cubic and hexag- 
onal sitesi^i A relevant effect for the diffusion in n-type 
4H-SiC is the stabilization of additional negative charge 
states (Vq, V 2 ^, V s j~, and V s j~) due to the larger band 
gap. 

For the migration of the vacancies we have analyzed 
mechanisms that involve either nearest neighbor hops or 
second neighbor hops. In Fig.[3|both mechanisms are in- 
dicated for the carbon vacancy. Our analysis has shown 23 
that the silicon and carbon vacancies migrate by second 
neighbor hopsiiS* 4 ^ The arguments that rule out a migra- 
tion entirely based on nearest neighbor hops are as fol- 
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lows. By a nearest neighbor hop the vacancy transforms 
into a vacancy-antisite-complex as the silicon (carbon) 
neighbor moves onto the sublattice of the opposite kind, 
thereby creating an antisite next to the vacant lattice 
site. In the case of the carbon vacancy we have found 
that the silicon vacancy-antisite-complex is unstable in 
all relevant charge states? 2 ^ A barrier that could stabilize 
the vacancy-antisite complex is not present. This find- 
ing rules out any migration mechanism based on nearest 
neighbor hops for the carbon vacancy. In case of the sil- 
icon vacancy, it has been found by us2£ in 3C-SiC and 
by Rauls et alm^ for the neutral complex in 4H-SiC that 
the carbon vacancy-antisite-complex may be more sta- 
ble than the vacancy. This metastability of the silicon 
vacancy occurs in p-type and intrinsic (compensated) 
material and has implications on the migration as well 
as the annealing of Vsi- It will be considered in detail 
in Sec. IIII CI However, a migration mechanism entirely 
based on nearest neighbor hops is impossible: a consec- 
utive nearest neighbor hop would transform the Vc-Csi- 
complex into a Vsi-Sic-Csi-complex. We have found that 
this complex is unstable. 

Our analysis of the second neighbor hop, its migration 
path and the calculation of the migration barriers is dis- 
cussed in detail in Ref. Here we basically summarize 
our findings for the migration barriers as tabulated in 
Tab IIVI The barriers are obtained using the 216 atom 
cell and the T-point for the k-summation. This also al- 
lowed for the proper inclusion of Jahn- Teller distortions 
in case of the carbon vacancy and the inclusion of spin 
polarization for the silicon vacancy. Using a 64 atom cell 
and special k-points (2x2x2 Monkhorst-Pack mesh) 
the calculated barriers agree within 0.3 eV. For the car- 
bon vacancy we find relatively high migration barriers. 
They strongly depend on the charge state and vary be- 
tween 5 eV in p-type material (Vq + ) and 3.5 eV in n-type 
material (Vq). The charge state dependence has its ori- 
gin in the successive occupation of bonding defect levels 
at the transition state. For the silicon vacancy we obtain 
smaller migration barriers and a less pronounced depen- 
dence on the charge state. The migration barrier in p- 
type and intrinsic material with values between 3.2 eV 
and 3.6 eV is larger than the barrier for the transfor- 
mation into the carbon vacancy-antisite-complex. The 
implications of this finding are discussed in the next sec- 
tion. In n-type material the migration barrier amounts 
to 2.4 eV (Vg~). Since the vacancy is stable in n-type 
material, a possible transformation of the vacancy here 
does not interfere with its migration. 



C. Metastability of the silicon vacancy 

Our results in 3C- and 4H-SiC for the metastability of 
the silicon vacancy (in 4H-SiC we have considered com- 
plexes with the vacancy-antisite pair located at neighbor- 
ing cubic or hexagonal sites and aligned along the c-axis) 
and the findings of Lingneref alm& suggest that the follow- 




FIG. 3: Migration of the carbon vacancy: (1) nearest neigh- 
bor hop (2) migration on the carbon sublattice by second 
neighbor hops. 



TABLE IV: Vacancy migration and transformation of Vsi into 
Vc-Csi: Energy barriers for the relevant charge states. 



path 




enerj 


ry barriers 


[eV] 






2+ 


1+ 





1" 


2" 


V c -*Vc 


5.2 


4.1 


3.5 






V S i^V S i 




3.6 


3.4 


3.2 


2.4 


Vsi^Vc-Csi 




1.9 


2.4 (2.2) 


2.5 


2.7 


Vc-Csi^Vsi 


6.1 


4.2 


3.5 (2.4) 


2.4 





ing discussion is not restricted to a certain polytype. The 
silicon vacancy and the carbon vacancy-antisite-complex 
may be transformed into each another by a nearest neigh- 
bor hop as discussed above. Depending on the doping 
conditions a metastability of the silicon vacancy arises, 
when the Vcj-Csi-complex is more stable than the silicon 
vacancy. This is the case in p-type and intrinsic (com- 
pensated) material. A stable silicon vacancy is retained 
for a Fermi-level position above 1.7 eV. This is demon- 
strated in Fig. 01 For a given Fermi- level, the two defects 
Vsi and Vc-Csi possess different charge states, as the 
transformation of Vsi into Vc-Csi rises the position of 
localized states within the band gap. 

The energy barriers we have obtained for the trans- 
formation Vsi — ►(Vc-CsO an d the reverse process (Vc- 
CsO^Vsi are listed in Tab. IIVI The values refer to our 
recent results obtained for 3C-SiC using a 216 atom cell 
and the T-point and include spin effects^ The barri- 
ers are calculated for each given charge state. We as- 
sume that the initial charge state is preserved during 
the transformation and equilibrates at the final con- 
figuration. This applies also to the transformation of 
Vf- ^(V c -Csi) 2 - and (V c -C S i) 2+ ^V 2 +, where the 
final configurations (Vc-Csi)~ and Vg ; are the ground 
state configurations. For all charge states the transfor- 
mation occurs in a low-spin state, except for Vg; where 
the high-spin state has a lower barrier (given in braces 
in Tab. If VI) . Since this high-spin state is not the ground 
state of the final Vc-Csi-complcx, the complex relaxes 
into the low-spin state before a reverse transformation 
occurs. 

The metastability of the vacancy in p-type and its 



7 



> 




R F [eV] 

FIG. 4: Formation energy of the silicon vacancy Vsi and the 
carbon vacancy-antisite-complex Vc-Csi for 3C-SiC and 4H- 
SiC. The experimental band edges of 3C- and 4H-SiC are in- 
dicated by vertical bars. The subscripts refer to the polytype 
and the cubic and hexagonal lattice site (k and h respectively) . 
(Vc-Csi)4H,fc and (Vc-Csi)4H,h refer to complexes with the 
vacancy and antisite at either cubic or hexagonal sites that 
are aligned along the c-axis. 



stabilization in n-type material is accompanied with a 
pronounced variation of the transformation barriers with 
the charge state. In p-type material the transformation 
Vsi — ► Vc-Csi has a much lower barrier with values be- 
tween 1.9 and 2.5 eV than the reverse transformation 
with a barrier between 6.1eV and 3.5 eV. The transfor- 
mation is thus more likely than the reverse transforma- 
tion. Only when (Vc-Csi) - becomes stable, i.e. inn-type 
3C-SiC, the transformation of V|j~ has a lower probabil- 
ity than the transformation of (Vc-Csi) - . 

A comparison of the transformation barrier with the 
migration barrier shows that the transformation occurs 
with a much higher probability than the migration in the 
charge states 1 + , neutral and 1 — , which are relevant in 
p-type and compensated material. The reverse transfor- 
mation Vc-Csi^Vsi has much larger barriers than the 
migration of Vsi- Note that the relevant charge states are 
2 + and 1 + for this process in p-type material. Therefore, 
the vacancy migration is kinetically hindered for low and 
moderate temperatures unless the reverse transformation 
is activated at high temperatures. Then a migration of 
Vc-Csi may take place via Vsi- The effective barrier 
for the migration path Vc-Csi— *Vsi -^Vsi — ►(Vc-Csi) 
is given by the formation energies E f Vsi , its migration 
barrier E^ 3i , and the formation energy 

E Jo-C si of the 

complex 



E, 



v Si 



E, 



Vc-C Si 



E v Sl 



(4) 



In p-type material this effective barrier amounts to 7.7 eV 
and drops to a value of 4.3 eV for /ip = 1.2 eV around 
mid-gap, when the vacancy is negative. The effective 
barrier is still 3.2 eV at — 1.76 eV, when the silicon 
vacancy becomes doubly negative. Here we have assumed 



TABLE V: Recombination of vacancies and interstitials: En- 
ergy barriers of Frenkel pairs for the relevant charge states. 
The geometry of the Frenkel pairs is displayed in FigQ>] 



Frenkel pair 



energy barriers [eV] 



C. 



p(100) 



C. 



spSi(lOO) 



Si 



sp(HO) 



-V c 
-Vc 

•Vsi 



2+ 
1.0 
1.2 

0.0 



V 







0.5 
1.2 



0.2 



0.4 
1.4 

0.2 



that the charge state equilibrates at the vacancy. If this 
is not the case, the effective migration barrier for the pos- 
itive and neutral Vc-Csi-complexes are 6.3 eV and 4.5 eV 
respectively. 

Besides a migration of a metastable silicon va- 
cancy, also a dissociation of the carbon vacancy-antisite- 
complex has to be considered^ In this case the carbon 
vacancy migrates away from the antisite at the expense 
of the binding energy. Using the migration barrier for 
V c between 3.5 eV (V c ) to 5eV (V c + ) plus the binding 
energy of about 1 eV as an estimate for the dissociation 
barrier, we obtain values of 6.2 eV for (Vc-Csi) 2+ and 
4.5 eV for (Vc-Csi)°. Only in p-type material this pro- 
cess has a similar probability as the reverse transforma- 
tion followed by a migration of the silicon vacancy. 

We have not investigated the transformation barriers 
in 4H-SiC. Rauls et alm& have obtained a transformation 
barrier of 1.7 eV for the neutral vacancy using a DFT- 
bascd tight-binding scheme, which is in agreement with 
our results. Therefore we expect a similar charge state 
dependence of the barriers in 4H-SiC as discussed above. 
For the charge states 3~ and 4~ of the vacancy we ex- 
pect larger barriers for the transformation Vs; — ► (Vc- 
Csi) and correspondingly lower barriers for the reverse 
transformation. 



IV. RECOMBINATION OF FRENKEL PAIRS 
A. Carbon vacancies 

In a carbon Frenkel pair the interstitial and vacancy 
may be nearest neighbors or second neighbors. Also 
vacancies and interstitials may attract each other over 
larger distances. For their recombination the interstitial 
has to migrate to the vacancy. In the following we ana- 
lyze the recombination of Frenkel pairs where a carbon 
split-interstitial is a nearest or second neighbor of the 
carbon vacancy. 

The nearest-neighbor-pair C sp si(ioo)~Vc is a vacancy- 
complex with a carbon split-interstitial at the silicon site. 
Two different configurations are possible: a stable con- 
figuration with a threefold coordinated carbon intersti- 
tial and a silicon dangling bond as depicted in Fig [SJi, 
and another configuration in which the carbon interstitial 
possesses one carbon neighbor. This latter pair turned 
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(a) 



(b) 



(c) 




FIG. 5: Recombination of carbon and silicon Frenkel pairs: 
(a) carbon nearest-neighbor-pair C sp si(ioo)-Vc and second- 
neighbor pair C sp (ioo)-Vc, (b) silicon Frenkel pairs with Sitc- 
interstitial (the labels SiTSi and SItc, i indicate the unstable 
sites) and (c) fourth-neighbor-pair Si sp (iio-Vsi- The lowest 
energetic path of the interstitial atom is indicated by arrows. 



out to be unstable. The first configuration exists in the 
charge states 2 + through 1~ with the ionization levels 
at 0.4eV (+|2+), 0.7eV (0|+), and 1.8 eV (-|0). The 
defect states are formed by a silicon p-orbital of the in- 
terstitial and the dangling bond of the vacancy within the 
same plane. Using the saddle point search method and 
a 64 atom cell with special k-points we have obtained 
a recombination barrier for the neutral charge state of 
1.4 eV (c.f. Tab. Iv|) . At the transition state the occu- 
pied defect levels are well localized. Towards the end 
of the recombination they transform into dclocalized va- 
lence band states. The pair may also recombine in the 
positive charge states. There we expect the neutraliza- 
tion by an electron transfer at the end of the recombi- 
nation. In this case a recombination barrier can be cal- 
culated and we obtain 1.2 eV for (C spSi { 10 o)-Vc) 2+ and 
(C spS i(ioo}-Vc) 1+ - For (C spSi(10 o)-V c )~, which is rele- 
vant in n-type material, a second defect level is occupied 



that transforms into a delocalized conduction band state 
during the recombination. Two possibilities exist in this 
case: firstly, the defect is converted to neutral by an elec- 
tron transfer before the recombination or secondly, in the 
final stage of the recombination the additional electron 
is donated into the conduction band. The treatment of 
both cases is beyond the scope of the present work. 

In the second-neighbor-pair C sp (ioo)-Vc the interstitial 
is located at the neighboring carbon site, as depicted in 
Fig. The pair exists in the charge states 2 + , 1 + , 
and l - , with the ionization levels e(+|2+)=0.8cV, 
e(0|+)=l.leV and e(-|0)=1.6eV. The defect states 
mainly derive from the interstitials levels. The recom- 
bination occurs by a second-neighbor hop of the upper 
carbon atom via the unstable C sp si(ioo> -configuration as 
indicated in Fig. [Hp,. The recombination barriers for the 
pair are listed in Tab.lVl With values between 1.0 eV and 
0.4 eV they are considerably lower than the barriers for 
the recombination of C sp gj{ 10 o)-Vc. The barrier is also 
lower than the migration barrier of the split-interstitials. 
For the negative charge state a similar mechanism as for 
the negative nearest-neighbor pair has to be considered. 

In p-type material the vacancy and the approaching 
interstitial may be positively charged. An additional long 
range Coulomb repulsion may hinder the recombination. 
The barrier for the formation of the Frenkel pair is then 
given by the migration barrier of the interstitial plus the 
Coulomb repulsion of the pair. On the other hand, once 
the defect levels begin to overlap new states are formed. 
The charge state will be successively reduced by electron 
transfer to the defect levels. Indeed, for the Frenkel pairs 
discussed above the total charge of the isolated defects 
Cwioo) and is reduced to 2+ or 1 + . 

A considerable energy gain is associated with the for- 
mation of the Frenkel pairs C sp gi/ioo)-Vc and C sp ^ 10 o)- 
Vc from the isolated defects C sp ( 10 o) and Vc- The 
binding energy for C sp si(ioo)-Vc varies between 2eV 
for (C spSi<100) -V c ) 2 + ^C° p(100) +V 2 + and 4eV for 
(C spS i<ioo>-Vc)° ^C s ° p(100) +V£. For the second- 
neighbor pair C sp (ioo)-Vc we obtain values between 
0.8 eV for (C sp(100> -V c ) 2 + -C° p(100) +V 2+ and 1.0 eV 
for the separation of the neutral pair. The interstitial 
is thus attracted by the vacancy once the charge state 
has equilibrated. 

Recently, the recombination of vacancies with car- 
bon split-interstitials was investigated in 4H-SiC using a 
DFT-based tight-binding schemed Only neutral Frenkel 
pairs were considered. A similar recombination barrier of 
the neutral pair C sp ( 10 o)-Vc (0.5 eV) was found as in our 
calculation for 3C-SiC. Larger barriers were obtained for 
the direct recombination of more remote pairs. 



B. Silicon vacancies 

Frenkel pairs involving the silicon vacancy may be 
formed with either of the dominant interstitials Sixco 
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Si sp (iio) and with the silicon-coordinated interstitial 
Sixsi- The Frenkel pairs with Sixc an d Sixsi should ac- 
count for the recombination in p-type material, when the 
SiTC-site is relevant for the migration. A recombination 
via the Frenkel pair with Si sp ^ 110 ) should be relevant in 
intrinsic (compensated) and n-type material, where split- 
interstitials dominate. 

First we have analyzed Frenkel pairs formed with a 
SiTSi- or SiTC-interstitial neighboring the silicon vacancy. 
The Frenkel pairs with a Sixc-interstitial or a Sixsr 
interstitial next to the vacancy turned out to be unstable 
in all relevant charge states (the labels Sixc, 1 and Sixsi 
indicate the sites of the interstitial in Fig. |SJd). Upon 
reaching these sites, the interstitial immediately recom- 
bines with the vacancy. The closest stable Frenkel pair 
we found involves an interstitial at a Sixc-she with only 
one common carbon neighbor to the vacancy (c.f. the 
Sixc, 2-site in Fig.Eb). This stable Frenkel pair has deep 
levels within the band gap. It is only realized in the 
positive and neutral charge state with ionization levels 
at e(+|2+) = 0.5 eV and e(0\+) = 1.2 eV. The recom- 
bination of this pair could be achieved via a hop of the 
interstitial passing close by the adjacent Sixsrsfte. Our 
calculations indicate that this is not a likely path. In- 
stead an anti-structure pair forms by a hop of the carbon 
neighbor into the silicon vacancy. In a concerted mo- 
tion the interstitial moves into the site left by the carbon 
neighbor. The reaction barrier amounts to 3.2 eV for the 
positive and the neutral Frenkel pair, which is compara- 
ble with the migration barrier of Sixc • 

Frenkel pairs with the split-interstitial Si sp (no) rep- 
resent other possible configurations. It turns out that 
Si sp (no) -Vsi is unstable when the vacancy and the in- 
terstitial are second neighbors, irrespective of the orien- 
tation of the pair. A recombination barrier thus may 
be encountered for more distant pairs, e.g. with the 
interstitial at the third or fourth neighbor shell. The 
recombination of such pairs is based on the migration 
of the split interstitial to an unstable second neighbor 
position. Here we discuss the fourth-neighbor Frenkel 
pair in more detail. Our findings for this pair suggest 
a similar mechanism for the third neighbor pair. The 
fourth-neighbor complex is stable except in highly doped 
p-type material. In the doubly positive charge state the 
interstitial and the vacancy immediately recombine. The 
complex becomes positively charged at the ionization 
levels (2 + |l + )=0.48eV. For a Fermi level above mid- 
gap the complex is neutral or negative (ionization lev- 
els: (l+|0)=1.3eV and (0|l-)=1.8eV). It is doubly neg- 
ative in strongly doped n-type material only (ionization 
level(2 _ |l~)=2.3eV). The saddle point search finds the 
following recombination mechanism (c.f. Fig. |Sp): first 
the silicon interstitial jumps towards the silicon-neighbor 
it has in common with the vacancy, then this neighbor 
is kicked out of its site and recombines with the vacancy. 
This process is accompanied with a barrier of 0.2 eV. As 
the split-interstitial close to the vacancy is unstable, the 
barrier for this process is considerably lower than the 



TABLE VI: Dissociation energy of neutral interstitial and 
vacancy clusters. For the interstitial clusters the given value 
is the energy needed to remove a single carbon atom. 



Poly type 


Cluster 


Dissociation energy (eV) 


3C 


( W«n ) 9 


2.8 


3G 


(C sp ) 3 


3.0 


3C 


(C sp )4 


5.7 


4H 


(C sp )2,kh 


5.8 


4H 


(C S p)2,hh 


5.0 


4H 


(C S p)2,kk 


4.6 


3C 


Vc-Vsi 


4.5 


3C 


Vc-Vc 


1.7 


3C 


Vsi-Vsi 


0.1 



TABLE VII: Ionization levels of the investigated interstitial 
clusters. The tri-interstitial (C sp )3 exhibits a negative-?/ ef- 
fect for the positive charge states and the tetra-interstitial 
(C sp )4 is electrically inactive. 







£ (1+|2+) 


£ (0|1+) 


£ (1-|0) 


e(2-\l-) 


3C 


(Csp)2 


0.8 


0.99 






3C 


(C sp ) 3 


0.85 


0.62 


1.49 


1.8 


3C 


(C S p)4 










4H 


(C S p)2,kh 




0.10 


2.78 


3.13 


4H 


(Csp)2,lih 


0.41 


0.50 


2.82 


2.68 


4H 


(C sp )2,kk 


0.32 


0.38 


2.66 


2.6 



barrier for the interstitial migration. As discussed for the 
annealing of carbon vacancies, we expect a higher barrier 
for the recombination in the negative charge state of the 
complex. 

Thus the analysis of the silicon Frenkel pairs shows 
that the pairs with a short distance between the inter- 
stitial and the vacancy are unstable, except for the sec- 
ond neighbor pair pair Vsi-Sixc.2- For the stable pairs 
we have found recombination barriers that are some- 
what lower than the migration barriers of the interstitial. 
These findings suggest a recombination based on a migra- 
tion towards the closest stable sites. They also indicate 
that the energy barriers of this process should not exceed 
the migration barrier of the corresponding interstitial. 



V. CLUSTERING OF CARBON 
INTERSTITIALS 

A formation of composite defects, e.g. Di or Dn cen- 
ters, during annealing requires a reservoir of intrinsic de- 
fects. However, at elevated temperatures vacancies and 
interstitials anneal out. (see discussion in Sec. I VIII . Yet, 
the point defects can be kinetically supplied by defect 
precipitations. The clusters of vacancies and interstitials 
can then serve as sources of mobile defects. We expect 
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FIG. 7: Di-interstitial (C sp )2.kh in 4H-SiC located on a cubic 
and a hexagonal site. The subscripts k and h refer to cubic 
and hexagonal sites. Ci refers to the additional carbon atoms. 
The stacking sequence of the crystal is indicated by the dark 
bonds. 



that the effect of the interstitial clusters on the defect 
kinetics is much more pronounced than that of the va- 
cancy clusters. Although both cluster types possess size- 
able dissociation energies (cf. Tab. IVII and Ref. |^ for 
the divacancy), the migration barriers of the carbon and 
the silicon vacancies are much higher than those of the 
interstitials (cf. Sec. IIIII) so that it is much more likely 
that the interstitials combine into precipitates. In addi- 
tion, we expect (cf. Sec. IVIfl that the vacancy- interstitial 
recombination sets in before a significant vacancy clus- 
tering occurs. 

We focus in this section on small clusters formed by 
upto four carbon interstitials. Due to the high mobility 
of the carbon interstitials, these clusters may act as a 
sink for interstitials at lower temperatures and can re- 
emit them at higher temperatures. Similar aggregates 
have already been investigated in diamond^ The micro- 
scopic structure of the considered complexes is displayed 
in Fig. H3 for 3C-SiC, an example of the di-interstitial in 
4H-SiC is shown in Fig. [7| A detailed analysis of the 
electronic structure and the vibrational spectra is given 
elsewhere Here we summarize the facts that are im- 



portant for the annealing mechanisms. The focus is on 
the dissociation energies of these defects, i.e. the energy 
needed to remove a single carbon atom from these clus- 
ters. For a mid-gap Fermi-level all clusters are neutral. 
The results for the dissociation energy in the neutral 
charge state are listed in Tab. IVII 

The carbon split-interstitial possesses two dangling p- 
orbitals oriented in (llO)-direction, resulting from the 
sp 2 -hybridization of the carbon atoms. The two adja- 
cent split-interstitials form a covalent bond due to the 
overlap of these orbitals. This (cf. Fig. HJi) results in 
a more favorable sp 3 -hybridization of the two interstitial 
carbon atoms and gives rise to an energy gain of 2.8 eV for 
the carbon di-interstitial (C sp )2- In the same way further 
carbon atoms can be absorbed, resulting in larger carbon 
aggregates. The tri-interstitial (Fig. EJd) possesses a sim- 
ilar dissoziation energy (3.0 eV) as the di-interstitial. In 
the tetra-interstitial a ring of split-interstitials is formed 
(Fig. which is especially favorable. All carbon atoms 
are sp 3 -hybrids, hence all bonds are saturated. This leads 
to an energy gain of 5.7 eV for the added carbon intersti- 
tial. 

The charge states of the investigated clusters are listed 
in Tab. I VI il The di-interstitial possesses all charge states 
from 2+ to 0. The tri-interstitial is mainly neutral, only 
at the band edges the charge states 2+ and 2— can exist. 
This cluster exhibits a negative-?/ effect, since it is ener- 
getically favorable to fully occupy the connecting bonds 
between all three adjacent split-interstitials (cf. Fig.[U>). 
The tetra-interstitial is electrically inactive, since all car- 
bon p-orbitals are saturated. 

In 4H-SiC, a stronger relaxation and rehybidization 
of the carbon atoms is observed than in 3C-SiC. There- 
fore already the di-interstitial shows a very high stability. 
Three different variations of the di-interstitial are possi- 
ble. The split-interstitials can either lie completely in a 
hexagonal or a cubic plane, or they occupy two adjacent 
cubic and hexagonal sites (cf. Fig. 01. With a binding 
energy of 5.8 eV the latter is the most stable configura- 
tion. This higher binding energy compared to 3C-SiC 
results from the sp 3 -hybridization of the border atoms, 
which is preferred due to the different geometry of the 
hexagonal planes in 4H-SiC. Also the hexagonal and cu- 
bic di-interstitials are very stable with binding energies 
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of 5.0 eV and 4.6 eV, respectively. These di-interstitials 
show a strong distortion and have a nearly linear struc- 
ture within their cubic or hexagonal plane. 

The micro-clusters presented in this section are exam- 
ples of clustering processes. The formation of cluster net- 
works via bonds between the p-orbitals of adjacent (C sp )2 
and (C sp )3 clusters is also possible. Additionally, carbon 
antisites can trap carbon interstitials with sizable binding 
energies 4& In the discussion of the annealing mechanisms 
below we outline the competition of the cluster formation 
with the vacancy-interstitial recombination and their role 
as a soure for carbon interstitials at the high temperature 
annealing stages. 



VI. HIERARCHY OF ANNEALING 
MECHANISMS 

The annealing of vacancies and interstitials is governed 
by several competing mechanisms. The principal mech- 
anisms are (i) the recombination of vacancies and inter- 
stitials, (ii) the out-diffusion to the surface, (hi) the dif- 
fusion to sinks, where they are annihilated or otherwise 
form stable complexes, and (iv) the transformation into 
other more stable configurations that is not diffusion- 
related. Sinks for interstitials or vacancies could be im- 
mobile defects, such as interstitial clusters and stable 
complexes with impurities, or extended defects, such as 
dislocations. In the latter cases the interstitial or va- 
cancy are not disappearing as such, but are bound in a 
more stable form. In the experiment, nevertheless, the 
observed signatures anneal out, while other centers re- 
lated to the annealing may appear in the spectra. 

The experimental investigation of the annealing kinet- 
ics requires a sufficient temperature T to activate one or 
more of the above annealing mechanisms. The diffusion- 
limited mechanisms (ii)-(iii) and the transformation- 
related mechanisms are typically described by first order 
kinetics with an activation energy given by the migra- 
tion or transformation barrier (e.g. Ref. 54). The ki- 
netics of the vacancy-interstitial recombination depends 
on the average separation of the vacancy and intersti- 
tial. For bound Frenkel pairs the activation energy is 
given by the recombination barrier. For uncorrelated 
vacancy-interstitial pairs (e.g. when their separation is 
much larger than the capture radius) the interstitial has 
to migrate to the vacancy (the vacancies are rather im- 
mobile) to form a Frenkel pair. The Frenkel pair then 
simply recombines by a hop of the interstitial into the 
vacant lattice site. In this the recombination bar- 

riers of the Frenkel pair are found to be relatively small, 
the recombination is essentially limited by the diffusion of 
the interstitial. The recombination then is described by 
the activation energy of the interstitial migration. The 
prefactors of the diffusion-limited mechanisms and the 
vacancy-interstitial recombination, however, depend on 
parameters arising from the sample prepartion, e.g. the 
average distribution (and concentration) of the defects 



(and sinks), as well as the mutual interaction between 
the involved partners. Hence, the predictive translation 
of activation energies into annealing temperatures hinges 
on parameters that are specific to the particular sam- 
ples. Only by detailed simulations it is possible to assess 
the initial state of the samples (e.g. after irradiation) 
and the annealing kinetics quantitatively. Our analysis 
forms the basis for such simulations, that are well be- 
yond the scope of the present paper. Nevertheless, a 
qualitative assessment of the kinetics via the calculated 
activation energies should be valid, since the activation 
energies differ substantially in value and enter the rate 
constants exponentially. 

In the following we deduce a hierarchy of annealing 
mechanisms based on a comparison of the energy barriers 
involved in the principal annealing mechanisms. Central 
to the discussion are four important observations. Our 
results for the migration of intrinsic point defects show 
that the interstitials are by far more mobile than the va- 
cancies, also the recombination barriers for Frenkel pairs 
are lower than the migration barriers of the vacancies. 
The silicon vacancy is a metastable defect in p-type and 
compensated material. The migration of defects strongly 
depends on the charge state, in particular the migration 
of silicon interstitials and vacancies. 



A. Silicon vacancies and interstitials 

Qualitatively, the annealing of silicon vacancies is sub- 
ject to three different Fermi-level dependent mechanisms: 
(i) the transformation of the silicon vacancy into a carbon 
vacancy-antisite complex, (ii) the Frenkel pair recombi- 
nation with silicon split-interstitials and (iii) the migra- 
tion of silicon vacancies to sinks. These mechanisms are 
summarized in Fig. ^ 

In p-type material, we expect a dominance of the 
annealing of the vacancy signature by its transforma- 
tion into the vacancy-antisite complex. The recom- 
bination of the stable Frenkel pairs, i.e. pairs with 
the Sixc-interstitial located in the second neighbor 
shell or beyond, is suppressed by the larger migra- 
tion/recombination barrier (3.2eV-3.5eV) as compared 
to the barrier of the transformation (1.9eV-2.4eV). Also 
the migration of the silicon vacancy is kinctically sup- 
pressed in p-type material due to the larger migration 
barrier (3.6eV-3.2eV). Once the silicon vacancy is trans- 
formed into the vacancy-antisite complex, it stays in this 
configuration. The reverse transformation requires fairly 
high temperatures due to the large barrier of 6.5 eV. As 
discussed in Sec. IIII CI the dissociation of the vacancy- 
antisite complex becomes possible at the same time. 

In compensated material (with the Fermi level at or 
above mid-gap), the silicon split-interstitial is the most 
relevant interstitial configuration. Its migration barrier 
is lower than the transformation barrier of the silicon 
vacancy. Hence the vacancy-interstitial recombination 
prevails in the annealing hierarchy. For the recombina- 
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tion of bound Frenkel pairs we found a low barrier of 
0.2 eV. The annihilation of vacancies by more distant in- 
terstitials is limited by the migration of the interstitial 
towards the vacancy with a migration barrier of 1.4 eV. 
It constitutes a separate annealing stage. Note that the 
small binding energy of the Frenkel pair indicates a weak 
attractive interaction between the interstitial and the va- 
cancy. Hence a dissociation of distant pairs has a rela- 
tively high probability, which enables further annealing 
stages. The metastability-induced annealing is the next 
annealing stage in the hierarchy. As a final stage, the 
vacancy-antisite complex, that evolves from the silicon 
vacancy, anneals by dissociation or by its migration to 
sinks. Again the migration of silicon vacancies does not 
directly contribute as the migration barrier of 3.2 eV to 
3.4 eV (Vgj and Vg 1 ; respectively) is larger than the trans- 
formation barrier. 

In n-type material the silicon vacancy becomes a stable 
defect. As in intrinsic (compensated) material the an- 
nealing by vacancy-interstitial recombination has a lower 
activation barriers than the vacancy migration. The mi- 
gration of silicon vacancies to sinks is most likely the last 
stage in the hierarchy. It has a slightly higher barrier 
than the transformation Vs; — >Vc-Csi in compensated 
material. 

By electron spin resonance experiments Itoh et alH 
observed the annealing of the Tl-center in p-type and n- 
type 3C-SiC irradiated by electrons and protons. They 
found three annealing stages at 150°C, 350°C and 750°C. 
The annealing behavior was found to be insensitive to 
the doping (aluminum and nitrogen) and the irradiating 
particles (electrons and protons), except for the p-type 
electron-irradiated sample, where the final annealing oc- 
curred at 350°C in coincidence with the annealing of 
three other defect-centers. Itoh et al. argued that a pos- 
sible loss of the paramagnetic state might be correlated 
with the annealing of two unidentified centers T6 and T7 
that were only observed under these conditions i 7 . The 
independence of the annealing behavior of the particular 
dopants indicates that the observed annealing process it- 
self involves only intrinsic defects. 

The Tl-center was identified by the authors as a 
negative silicon vacancy Vg;. Later this identifica- 
tion was verified theoreticallyi^iSLSiiSi by calculations of 
the hyperfine tensors. According to earlier theoretical 
resultsiasiiaSi and our calculations (c.f. Tab.[lJ the silicon 
vacancy is in a negative charge state over a wide range 
of the Fermi level starting below mid-gap (fiF >0.6eV) 
and ending somewhat below the conduction band edge 
(/if <1.76eV). This suggest that the p-type and n-type 
samples were compensated by the irradiation, most likely 
the Fermi level was trapped by irradiation-induced deep 
defects around mid gap. Under these conditions the 
scenario for compensated material we described above 
should explain the annealing of the EPR-signature of 
Vgj by three annealing stages, namely: (i) the initial 
recombination of Frenkel pairs with a small separation 
between the vacancy and interstitial, (ii) the recombina- 



tion of vacancies and interstials with a larger separtion 
that is limited by the diffusion of the interstitial, and (iii) 
the transformation of Vg; into a carbon vacancy-antisite 
complex. 

From the annealing experiments Itoh et a/S deduced 
an activation energy of 2.2 eV for the last annealing 
stage. According to our calculations the transformation 
of Vs; into Vc-Csi, that constitutes the third annealing 
stage, is associated with an activation energy of 2.5 eV 
(c.f. Tab. IIVI) in good agreement with the experimen- 
tal result. Once this transformation is activated silicon 
vacancies are largely transformed into carbon vacancy- 
antisite complexes and the EPR-signal of the Tl-center 
vanishes. During the short observation times (5 min) of 
the isochronal annealing the reappearance of Vg; is ki- 
netically hindered. Even though a signature of the car- 
bon vacancy-antisite complex may develop in the spectra, 
Itoh et al. reportedly have not observed the appearance 
of a new defect center that could evolve from the silicon 
vacancy. Note, that for a Fermi level below 1.24 eV the 
vacancy-antisite complex is not paramagnetic. Only in 
a small window, for a Fermi level between 1.24 eV and 
1.79 eV the paramagnetic positive charge state should 
prevail. Recently, Lingner et al. identified the Vc- 
Csi-complex by spin resonance experiments in neutron- 
irradiated n-type 6H-SiCiS& The hyperfine structure of 
an excited high spin state (S = 1) of this complex was 
detected using the magnetic dichroism of the adsorp- 
tion (MCDA) and MCDA-EPR. The theoretical analysis 
showed that the excited state is related to an intra-defect 
excitation of (Vc-Csi) 2+ , which is not paramagnetic in 
its ground state and therefore could not be detected by 
standard EPR. This is consistent with the observation of 
Itoh et al. The complex appeared only in samples an- 
nealed above the temperature at which the Tl-center in 
3C-SiC finally vanishes. EPR-centers with a similar g- 
tensor and fine structure constant D were observed ear- 
lier in n-type 6H-SiC^ 6 (P6- and P7-centers) and in elec- 
tron irradiated nominally undoped 3C-SiC (LS-center)^ 
The similarity suggests that the carbon vacancy-antisite 
complex as a common model for all these centers. Indeed, 
the annealing study of Son et al. shows that the center 
appears only after an annealing at 700°C-750°C. These 
observation support our interpretation of the annealing 
stages for the silicon vacancy reported by Itoh et al. . An 
diffusion-based annealing mechanism, that was proposed 
earlier^iS 7 . is apparently not involved in these experi- 
ments. 

Using PAS Kawasuso et alm&H have observed the an- 
nealing of vacancy-related defects in electron-irradiated 
n-type 3C-SiC and 6H-SiC. Two lifetime components 
were extracted from the positron signal in 6H-SiC. One 
component was attributed to defects related to silicon 
vacancies The other was shown to originate from 
carbon vacancy-related defects. In 3C-SiC only the com- 
ponent corresponding to silicon vacancy-related defects 
was detected. Kawasuso et al. found that this lifetime 
component annealed in several stages. For annealing 
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temperatures below 500°C two annealing stages were ob- 
served in 3C- and 6H-SiC that followed the findings for 
the Tl-center in 3C-SiC. In 3C-SiC vacancy-related de- 
fects were not detected above the annealing temperature 
of the Tl-center. This finding is also consistent with our 
interpretation of the annealing stages, since a positive 
carbon vacancy-antisite complex would repel the positron 
and therefore should not contribute to the positron life- 
time. 

In 6H-SiC, however, the signal persisted the anneal- 
ing at 750°C and finally annealed at 1450°C. Yet, an 
annealing stage with a less pronounced drop of the va- 
cancy concentration is present at 750° C. In a more re- 
cent PAS-study 6 using electron irradiated n-type 6H- 
SiC Kawasuso et al. describes the annealing of the 
silicon vacancy or related defect complexes with the a 
major stages at 500°C-700°C and 1000°C-1200°C. The 
vacancy-signature reached the bulk-level at 1200°C. In 
comparison to this more recent study, the samples of the 
earlier study had a larger nitrogen concentration (car- 
rier density of 5.5 x 10 17 cm -3 at room temperature vs. 
5 x 10 15 cm -3 ) and were irradiated with more energetic 
electrons at a lower dose (1 x 10 17 e~/cm 2 at 3MeV com- 
pared 3xl0 17 e~/cm 2 at 2MeV and). The differences in 
the annealing behavior is most likely related to these vari- 
ation in the nitrogen concentrations and the irradiation 
conditions. 

In the following we comment on presence of the anneal- 
ing stages above 1000°C in 6H-SiC. In their analysis of 
the earlier experiment Kawasuso et al. explained the final 
annealing stage at 1450°C in terms of the dissociation of 
vacancy-nitrogen complexes and the subsequent anneal- 
ing of the vacancy by out-diffusion or diffusion to sinks. 
The Vsi-N complexes were expected to form during the 
annealing. Such an interpretation is plausible as similar 
positron lifetimes are predicted for the Vsi-N complexes 
and the silicon vacancy^ Very recent theoretical inves- 
tigations show that these complexes are thermally stable 
in agreement with EPR-experimentsSS Yet, also the 
transformation into the carbon vacancy-antisite complex 
may play a relevant role in the PAS experiments, pro- 
vided that the sample remains compensated (Fermi level 
around 1.7 eV) and that the positron annihilation at the 
carbon vacancy-antisite complex resembles that for the 
silicon vacancy. Precise calculations based on ab ini- 
tio methods for this center are lacking sofar. However, 
calculations 59 using an approximate method based on the 
superposition of atomic densities (c.f. Ref. l6fj|) using the 
atomic coordinates obtained for Vc-Csi obtained in our 
work indicate for 4H-SiC, that the positron lifetime asso- 
ciated with this center lie in the range similar to the value 
of the silicon vacancy (at the cubic site between 162 ps 
to 175 ps and at the hexagonal site 187 ps to 193 ps de- 
pending on the charge state). Thus the positron lifetime 
does not significantly change when the negative silicon 
vacancy transforms into the neutral hexagonal vacancy- 
antisite complex. The migration of the carbon vacancy- 
antisite complex is associated with larger barriers than 



the transformation and hence is activated at higher tem- 
peratures. This scenario is a plausible explanation of the 
observed annealing of the Vsi-related signature in the 
PAS experiments, since the complex at the hexagonal 
site may become neutral during the annealing of compen- 
sating centers in n-type 4H-SiC. Obviously, Fermi level 
effects are relevant. The annealing characteristics should 
therefore depend on the polytype (due to the different 
band gap width). Based on our present understand- 
ing, we do not exclude the formation of vacancy-nitrogen 
complexes as an explanation of the final annealing stage. 
The interpretation of the last annealing stage in terms of 
vacancy-nitrogen complexes, however, requires that the 
dopant is present in similar or higher concentrations than 
the vacancies at the onset of the annealing stage. 



B. Carbon vacancies and interstitals 

For the carbon vacancy the annealing hierarchy starts 
with the vacancy-interstitial recombination at the first 
annealing stage and ends with the diffusion of vacan- 
cies to sinks or its out-diffusion (c.f. Fig. [TJ. The re- 
combination of Frenkel pairs containing the carbon split- 
interstitial as the nearest or second neighbor of the va- 
cancy proceeds with lower recombination barriers than 
the vacancy diffusion. Similarly the migration barrier for 
the carbon split-interstitial is lower than for the carbon 
vacancy. For the vacancy-interstitial recombination, the 
migration of the interstitial towards the vacancy is the 
bottle neck, since recombination of the nearest or second 
neighbor Frenkel pairs proceeds with a similar or lower 
barrier. The annealing of the vacancy based on its diffu- 
sion therefore should be activated at much higher tem- 
peratures than the vacancy-interstitial recombination. 

In irradiated material, where interstitials and vacan- 
cies are present in similar concentrations, we expect the 
vacancy-interstitial recombination to represent the first 
annealing stage. Our finding of a relevant binding energy 
for all charge states of the nearest and second neighbor 
Frenkel pair suggests an attractive interaction that traps 
the interstitial in the vicinity of the vacancy. Hence, we 
expect that a recombination of these Frenkel pairs has 
a higher probability than their dissociation. However, 
the diffusion-limited recombination of the isolated defects 
may be hindered by a Coulomb repulsion of the charged 
defects. This occurs only in p-type material for a Fermi 
level below 0.8 eV (cf. Tab.QJ. Above this value the car- 
bon split-interstitial is neutral. Even though the split- 
interstitial becomes negative in n-type 3C-SiC, the va- 
cancy should remain neutral. For n-type 4H-SiC, on the 
other hand, our calculations show that also the carbon 
vacancy becomes negative as a consequence of the larger 
band gap. The additional Coulomb barrier in p-type 3C- 
SiC and p-type or n-type 4H-SiC may leave a consider- 
able fraction of vacancies unannealed in this first anneal- 
ing stage. At the same time the carbon split-interstitials 
could anneal by diffusion to other defects where ther- 
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mally stable defects complexes or defect clusters could 
be formed. The remaining vacancies anneal in a second 
annealing stage based on the vacancy diffusion. 

So far the isolated carbon vacancy was not detected 
by EPR-experiments nor by PAS experiments in irradi- 
ated 3C-SiC. A spin resonance center T5 was observed 
together with the Tl-center (Vg ; ) by Itoh et al. that 
annealed at 150°C. Its hyperfme (HF) signature origi- 
nates from four silicon atoms and it was therefore iden- 
tified as a defect on the carbon sublattice. However, the 
original assignment to the positive carbon vacancy had 
been revised to a complex of two hydrogen atoms with a 
carbon vacancv^S 1 * 6 ^ given the /^-symmetry of the cen- 
ter that is incompatible with the theoretical findings^ 8 *^ 
for the carbon vacancy. Based on the calculations of 
HF-parameters for these models and the carbon split- 
interstitial it was recently shown that only the split- 
interstitial may explain the experimental findings for 
the T5-center. 23 i 24 i 63 i 64 Such an identification of the T5- 
center has to be verified experimentally, as the carbon 
HF-tensors of the (lOO)-oriented carbon dumbbell have 
sofar not been observed. Yet, this model would explain 
the low annealing temperature of the T5-center. Accord- 
ing to our analysis the positive carbon split-interstitial 
is indeed a highly mobile defect, with a migration bar- 
rier of only 0.9 eV. A comparison of the ionization levels 
of this defect and Vs; (cf. Tab. [1} shows that Cjp^ 100 ^ 

and Vgj simultaneous exist in a paramagnetic state for 
a Fermi-level around 0.7 eV in the p-type samples which 
is in agreement with the simultaneous observation of the 
Tl- and T5-centeriA£ At these doping conditions the car- 
bon vacancy prevails in the non-paramagnetic doubly 
positive charge state. This explains why the carbon va- 
cancy was not detected in EPR experiments or by PAS 
(the positron recombination at the positive vacancy is 
suppressed by the Coulomb repulsion). Therefore the 
mechanism (vacancy-interstitial recombination or other 
diffusion-limited mechanisms) behind the annealing of 
the T5-center cannot be deduced from the present ex- 
periments. 

A competing process to the vacancy-interstitial recom- 
bination is the formation of carbon interstitial clusters as 
discussed in Sec. Due to their high dissociation ener- 
gies such clusters would emit interstitials only at elevated 
annealing temperatures. These interstitials then take 
part in further defect reactions. For the di-interstitial 
and the tri-interstitial in 3C-SiC a dissociation energy of 
3eV plus an additional barrier, which may be approxi- 
mated by the migration barrier of the interstitial (0.5 eV 
for C° p ^ 100 ^), is needed to remove a carbon atom. This 
activation barrier is only slightly lower than the migra- 
tion barrier of the carbon vacancy (4.1 eV for V+ and 
3.5 eV for Vq). With a dissociation barrier of 5.7eV the 
tetra-interstitial is much more stable and should emit 
interstitials at higher temperatures than the di- and tri- 
interstitial. A comparison of the dissociation energies 
with the activation energies of the other processes sug- 
gests that new carbon interstitials should be supplied 



at temperatures above 1000°C. Indeed, the concentra- 
tion of e.g. the photo- luminescence centers Di 2 and Dn, 4 
which are presumably related to antisite o 65 ' 66 and car- 
bon aggregates^ 7 - rises significantly at annealing tem- 
peratures above 1000°CA^ For both formation processes 
the availability of free carbon atoms is an important pre- 
requisite. 65 ' 67 Defects similar to the carbon clusters de- 
scribed here are known from other semiconductors. Be- 
sides the carbon clusters in diamondS one can also men- 
tion here the well known {311}-interstitial defects in sil- 
icon. These defects emit silicon interstitials and thereby 
drive the transient enhanced diffusion of dopants during 
the annealing. 68 

In contrast to 3C-SiC, the signatures of the carbon 
vacancy and interstial were reported in 4H-SiC. In ir- 
radiated p-type 4H-SiC, recently two EPR-centers (Ell 
with S=l/2 and EI3 with S=l) were identified by Son 
et al. (Ref. |69J) that possess a similar g-tensor as the 
T5-center in 3C-SiC. The analysis of the HF-signature 
showed that the defects are located on the carbon sub- 
lattice. Theoretical calculations indicate that the ob- 
served silicon HF-tensors are consistent with the theo- 
retical values for the positive and neutral carbon split- 
interstitial, respectively4LSi§i Similar to the T5-center 
the Ell and EI3 center anneal around 200°C. In p-type 
4H-SiC irradiated at 400°C Son et a/J&S observed by 
EPR-experiments the final disappearance of the EI5- 
center at 500°C. The center was tentatively assigned to a 
positive carbon vacancy. This assignment was verified by 
calculations of HF-tensors ^2£2JL and it was shown that 
the center is located at the cubic site. The low annealing 
temperature of the EI5-center is consistent with our dis- 
cussion of the vacancy-interstitial recombination in 3C- 
SiC. The observation of only a single annealing stage for 
the carbon vacancy in the EPR-experiments is also con- 
sistent with PAS experiments 1 - performed in irradiated 
6H-SiC. There it was found that the corresponding life- 
time component annealed at 500° C. 

Recently, Konovalov et aZ.jSi Zvanut and Konovalov 
and Son et aW* have studied EPR-centers in semi- 
insula ting 4 H-SiC. They identified the centers ID1 
(Refs. l2~M70h and EI5 (Ref. EH). The latter was already 
observed in irradiated 4H-SiC. The center ID1 resembles 
the EI5 center with respect to its extensor and HF-values, 
which suggest a common identification as a carbon va- 
cancy at the cubic lattice site^ Both centers posses a 
higher thermal stability than the EI5-center in irradi- 
ated SiC. For example Konovalov et al. have found that 
the center persists up to 850°Ci2S The material used by 
Konovalov et al. and Zvanut and Konovalov was high pu- 
rity, as-grown (0001) 4H-SiC wafers grown by the seeded 
sublimation method. The high temperatures (> 2000 °C) 
applied during growth suggest that intrinsic defects are 
present in the material close to their equilibrium abun- 
dance. This means that carbon vacancies by far outnum- 
ber carbon interstitials, due to the much lower formation 
energy of the vacancy as compared to the interstitial. 
Also, the concentration of carbon interstitial clusters is 
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predicted to be well below the vacancy concentration. 
Therefore the excess concentration of carbon vacancies 
in a material grown at high temperatures may reach its 
equilibrium concentration at lower temperature only by 
a mechanism based on the vacancy diffusion. The large 
migration barrier of the carbon vacancy - 4.1 eV for Vj 
in 3C-SiC - explains the high thermal stability of carbon 
vacancies in as-grown semi-insulating samples. 

VII. SUMMARY AND CONCLUSION 

The microscopic picture of the annealing mechanisms 
of vacancies and interstitials has been developed from 
theoretical investigations based on ab initio methods 
within the framework of DFT. We analyzed in detail 
the ground state configurations of Frcnkel pairs formed 
by carbon (silicon) vacancies and interstitials and cal- 
culated the recombination barriers. We investigated the 
capture (and emission) of carbon interstitials by (from) 
carbon interstitial clusters. Various clusters involving up 
to four carbon interstitials were considered and the en- 
ergy needed to emit single carbon interstitials from these 
clusters was calculated. 

A hierarchy of annealing mechanisms has been de- 
rived that places the vacancy-interstitial recombination 
and the carbon clustering into context with diffusion 
based mechanisms and the mestability of the silicon va- 
cancy. The recombination barriers and dissociation en- 
ergies were compared with the barriers of interstitial and 
vacancy migration as well as the transformation barrier 
of the silicon vacancy. The relevance of the Fermi-level 
effect was demonstrated for the annealing hierarchy of 
the silicon vacancy and interstitial. The effect originates 
from two distinct configurations of the silicon interstitial 
in p-type and compensated (or n-type) material and the 
stabilization of the silicon vacancy in n-type material. In 
both cases a strong variation of the migration and trans- 
formation barriers with the Fermi- level is observed, which 
changes the hierarchical ordering of the annealing mech- 
anisms. While in compensated material the vacancy- 
interstitial recombination and the diffusion-based anneal- 
ing of silicon interstitials are activated at lower tempera- 
tures than the metastability-related annealing of the sil- 
icon vacancy this ordering is reversed in p-type mate- 



rial. In n-type material the metastability-based anneal- 
ing is unavailable and a diffusion-based mechanism is ob- 
served instead. Regarding the annealing of carbon inter- 
stitials and vacancies, the highly mobile carbon intersti- 
tials were shown to drive the vacancy-interstitial recom- 
bination and the formation of thermally stable carbon- 
interstitial clusters at the early stage of the annealing. 
The comparably low mobility of the carbon vacancy per- 
mits its diffusion-based annealing only at elevated tem- 
peratures. In this case, the hierarchical ordering was 
found to be independent of the Fermi-level effect. The 
emission of carbon interstitials from carbon interstitial- 
clusters is associated with activation energies that exceed 
the activation energies of the vacancy migration. Thus 
the clusters can provide carbon interstitials even at tem- 
peratures at which isolated carbon vacancies have van- 
ished. 

The annealing hierarchy was discussed in the light of 
annealing experiments on EPR- and PAS-ccntcrs, sup- 
ported by the recent microscopic identification of these 
centers through experimental and theoretical analysis. 
A consistent interpretation of the present experimental 
data for 3C-SiC is facilitated by the proposed anneal- 
ing hierarchy. The available experiments and theoreti- 
cal work for 4H- and 6H-SiC suggest that our qualita- 
tive conclusions should (with limitations) apply to other 
polytypes. We hope to stimulate further experimental 
investigation of the defects' annealing kinetics with our 
analysis. In particular, the annealing stages of the silicon 
vacancy above 1000°C and the confirmation of the ten- 
tative assignment of the EPR-centers T5, Ell and EI2 
to carbon interstitials call for further clarification. Also 
the role of defect aggregates, as exemplified by carbon 
interstitials, needs experimental substantiation. These 
aggregates should provide the missing link between the 
kinetics of the primary defects and the thermally stable 
PL-ccntcrs. 
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